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This paper investigates the use of Luenberger observer for sensorless power 
control of brushless double fed induction machine (BDFM) in wind energy 
conversion systems, the control strategy for flexible power flow control is 
developed by applying flux-oriented vector control (technique), In order to 
estimate the rotor speed, an adaptive algorithm based on Lyapunov stability 
theory will be design. Finally, the analyzed and simulation results in 
MATLAB/ Simulink platform confirmed the good dynamic performance of 
this new sensorless control for BDFG based variable speed wind turbines. 
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1. INTRODUCTION 

In recent years, the electrical machine has expanded considerably with the development of power 
electronics and data processing, in this way a many researchers developed the difference observation, for 
estimating the rotor speed and parameters identification of electrical machine. The Brushless double fed 
induction motor is one of the most important ac machines used because of its low cost and 
high reliability [1]. 

Sensorless control has been successfully applied to the BDFG based on an extended Kalman filter 
observer [2], The rotor speed estimator is designed by a phase locked loop ignoring the power winding 
resistance [3], and MRAS observer scheme based on the stator current of the control winding (CW) yessed 
the a phase locked loop (PLL) is proposed by [3]. 

The Luenberger observer is a well-known method for the sensorless control of cage induction 
machines, there are few reports related to the use of Luenberger observer for sensorless control of DFIG 
[4-6], when has been proved to be a good compromise between accuracy and complexity, and is able to work 
at wide speed range [7]. 

This paper discussed of a novel sensorless vector control of BDFG using Luenberger observer (LO), 
the error between the observed value and the true value considered the rotor speed , based on Lypunov's 
stability theory. General scheme of Luenberger observer is shown in Figure 1. 
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Figure 1. General scheme of Luenberger observer for speed estimation of BDFG 


2. THE MATHEMATICAL MODE OF BDFM 

The BDFM is normally operated in the synchronous mode and the natural synchronous 
speed equal to: 


Pp+Pc 


co p and co c are the angular frequency of power winding and control winding and Pp and Pc are the number of 
pole pairs of power winding and control winding. 

The mathematical model of BDFG with d-q reference (PW) is able to be expressed as [1], [8], [9]: 


Vp=Rp ' ip+ ^f +jC0p¥p 

V c = R c .i c + ^ + j(“p -( p p +PcK)V c 

V r =R r .i r +^ + j(co p -P p co r ) ¥r 

And the flux equations are given as: 

Vp=L p i p + M p i r 

Vc= L cic+ M cir 
Vr = L r i r +M c i c +M p i p 

The electromagnetic torque of BDFG is expressed as [3]: 

T e = | PpMp (i qp i dr - i dp i qr )-1P C M C (i qc i dr - i dc i qr ) 
The active and reactive powers of BDFM are as follows: 
Pp=f(V dp i dp +V qp i qp ) 

3 

Qp — 2 ^Vqp^dp Vlp^dp) 


( 2 ) 

(3) 

(4) 


(5) 

( 6 ) 

(7) 


( 8 ) 


(9) 

( 10 ) 
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3. VECTOR CONTROL DESIGN FOR BDFG 

In this section, the vector control of BDFM will be presented, to achieve regulation of the active and 
reactive power between the BDFG and the grid [9], [10]. The vector control of BDFM is similar of the 
principle of classical vector control of DFIM, which it based of annulled the quadrature component of the 
PW flux, and suppose the R p is neglected, the Equations (2) and (5) can be written as follow: 


v dp=° 


v qp = Vp=oy|/ P 


( 11 ) 


JVp=L p idp+ M pidr 

[0 = Lpi qp + M p i qr 


( 12 ) 


The rotor currents can be described using the power stator current: 

\J/p_ L^ 

L„ 


_ , A "''dp • 

ldr_ M“ _ M“ ldp 

p p 


M„ 


(13) 


3.1. The PW currents regulation 

The mathematical mode of BDFG in the steady state given by [11], 


[Vdp-Rpidp OJp. Lpiqp ®p Mpiqr 

[Vjp = Rp Jqp + COp.Lpijp + 0JpMpi c | r 


(14) 


Vic Rc-ldc ^pRclqc tOpM c iq r 

Si Si 

Vqc Rc’iqc ®pRcldc "L ^p^ddr 

S 1 S 1 


(15) 


0 = - Rr-idr - ® p L r i qr - co p M c i qc - o) p M p i qp 


0 = — R r i qr + ® p L r i dr + to p M c i dc + ® p M p i dp 


(16) 


Si, s 2 are the slips, which can be expressed as: 


m p-P|, m p c tOc-PptOp 


■, s 2 




OX 


(17) 


Used (14), (15) (16), (11) (13), The control winding can be expressed as 

V p L r 


i = ( L - Lp - Mp )i 

dc ( m p m c * 


dp 




MptOpMc M p M c © pSl 


(18) 


L r L n - M n R r \i/ 

i = (_LJE-L)i + 

QC V , , , , / A QD 


R r L p 


M p M c 


MpM c © p s 1 MpM c co p s, 


dp 


(19) 
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3.2. The CW currents regulation 

Used the Equations (3), (6), (18), (19), the CW voltage can be regulation by the CW currents as: 


V -R i +n M n di dc 

Me — Uc + ( B c t i - 


M c R r L p di qp 


L r dt co p L r s,M p dt 
“(Pp +PcK)(L c i qc +M c (-ijM qp )) 

M P i 


M c M p di dp 
L r dt 


( 20 ) 


M 2 di M R L di d M M di aD 

V _ T) : , /t c \ q c c r p ap p c qp 

ac -Rc-iqc +( L c- i 1 T - 


L r dt co L r SjM dt 


P r 1 P 

L, 


dt 


- («)p - (Pp + Pc K )(L c i d e + M c (-^ - -f i dp )) 


( 21 ) 


M p M p 


The third term: - (co p - (P p + P c )co r )(L c i dc + M c (^- - i dp )) Shows another cross and the general block 

control diagram is shown in Figure 2. 


M p M p 



Figure 2. Control scheme for BDFM 


4. THE LUENBERGER OBSERVER 

Using the six-order model of the Brushless doubly-fed induction machine in fixed stator d-q axis 
reference frame with PW current, CW current and rotor current components as state variables. 

The dynamic model of the BDFM is given in (d-q) reference frame that is used in FO for 
state observation, the model is given [12-18]: 


JX = AX + Bu 

[y = cx 

T 

u = [v pd V pq V cd 0 o] 

The state vector is 

T 

^ tpd ^pq kd kq hd hq ] 


( 22 ) 


The system matrix A, the input matrix B and the output matrix C are given as: 


1 

0 


0 

0 


0 

1 

0 

0 

0 

0 


, B - [mv(AL) ] ^ A = [_ inv(AL) AR ] 
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Where: 
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co p Mp 

0 
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0 

-co 2 M c 
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0 co 2 L c 

Rc 
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0 
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0 co 3 M c 

o 
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Where: 

<»2=<»p-(Pp+PcH 

«3 =Mp-PpM, 

The Luenberger observer which estimates the all stator currents will be designed using the BDFM model. 


X = AX+ BU + L(Y - CX) (23) 

L: is the observer gain matrix 

X = f A- Lcl X+ BU + LY (24) 


The Luenberger matrix gain L is chosen so that the poles of the characteristic matrix 
AL= A + LC to be stable. So, all eigenvalues of AL should have negative real parts. 

The poles can be placed by solving the differential equation, thus the matrix gain L can be 
calculated by the function (PLACE) Pole placement technique in MATLB. 

4.1. Estimation of the rotor speed 

The estimation error of the state variable giving by: 


E = (A-LC)E +AAX (25) 

AA is the error between the two matrices as being exclusively caused by the error between the real and the 
estimated speed 


AA = (A - A) 

Assuming that ^ cor = ®r — cor 


(26) 
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0 

-al 0 

a2 

0 

a3 

al 

0 -a2 

0 - 

-a3 

0 

0 

-a4 0 

a5 

0 

a6 

a4 

0 -a5 

0 - 

-a6 

0 

0 

a7 0 

-a8 

0 

a9 

- a7 0 a8 

0 - 

a9 
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The elements of matric are showed in the appendix. 

The speed observer can be constructed based on Lyapunov’s stability theory. 

Assuming that the Lyapunov function is defined as: 

V = E t E + — Acor 2 (27) 

kL 


Where: 


£ = 


l pd 1 pd 

ldpd lpd 

0 

0 

0 

0 


Where kL is a positive constant, Assuming that Acor = cor - cor 
The application of the general adaptation mechanism 


■= k J 


E AA X dt 


(28) 


Where: 


cor = k I e pd (-al i pq + a2 i cq + a3 i rq ) + e pq (al i pd - a2 i cd - a3 i rd ) dt (29) 


Which means: 


cor k j a2(e pd i C q e pq i cd ) + al(e p q i pd e pd i pq ) + a3(e pd i rq e pq i rd )dt (30) 

We can neglect the values of -al,a3 because its low values compared to a2, the adaptation 
mechanism becomes: 


cor = k|a2(e pd i cq -e pq i cd )dt 


(31) 


Where k is a positive constant. 

Usually the following proportional and integral adaptation mechanism, in order to improve the 
response of the rotor speed estimation. 

cor = Kp(e pd i cq - e cq i cd ) + KiJ (e pd i cq - i cd ) (32) 


5. WIND TURBINE MODEL 

In this work a horizontal axis wind turbine is used, which the mechanical power of the wind can be 
derived as: 
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P = |c p R 2 pv 3 (33) 

Where p =air density, R=radius of Blades, v =wind speed and C p =power coefficient which can be 
derived as: 


21 

C p (^,p) = 0.5176 (— -0.4p-5)e + 0.0068^ 

X: 


(34) 


Where: 


1 _ 1 0.035 

X + 0.08P p^ +1 

The power conversion coefficient defined as: 

V 


(35) 


(36) 


Where co t =the turbine rotor speed. 

The wind turbine is normally characterized between Cp and y for the given values of pitch angle 
(p°) is as illustrated in Figure 3. 
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Figure 3. Wind turbine generator C p - X characteristics 


5.1. Pitch angle controller design 

The advantage of pitch angle control is mor efficiency in low wind, small variation in the pitch 
angle can give strongly influenced by of the blade respect to the direction of the wind or to the plane 
of rotation. 

Used the wind velocity n, the reference rotor speed for extracting the MPPT is obtained by: 

_ ^op t v 


The gearboxes in a typical wind turbine increase the speed of the generator by the relation 


COm - G©t (38) 

The pitch controller is employed to regulation the rotor speed at the maximum used the rotor speed 
measure and the reference speed, which can find by the Equation (36). 

A simple proportional-integral (PI) controller is used to regulation, this regulator followed by 
limitation to fixing the angle to between the maximum and the minimum angle as shown in Figure 4. 
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Figure 4. Pitch angle controller 


Figure 5. Sensorless vector control of the BDFM 
using Fuenberger observer 


6. SIMULATION RESULTS 

The senseless control developed has been implemented in a MATLB 7.0 /simulation, The BDFM 
used in this simulation model is 3Y-3Y connected and its stator winding is 2-6 poles. The machine 
parameters presented by J. Poza [3] are used in this simulation as showed in table 1. 

To verify the state estimation performance extensive simulation tests were carried out to compare 
the sensouless control under different wind speed. 

A step change in wind speed is simulated in Figure 6, the wind speed is start at 5m/S, at 7second, 
the wind speed suddenly become 7m/S. 
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Figure 6. Wind speed 


Figure 7. Zoom of actual and estimated rotor speed 
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Figure 8. Error of rotor speed 
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Figure 9. Power winding reactive power 


Sensorless Control of Brushless Doubly-fed Generator Using Luenberger... (Hicham Serhoud) 

















































































































196 n 


ISSN: 2252-8792 



rr 

| 

-!-1 -!- 

:Lr| , 


i : i 


1 

| | | 

fTZ. 



1 i 

1 

.:.r. ; 

i i : 

*3 - a a w 


Figure 10. Power winding active power 
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Figure 11. Phase power winding current 




Figure 12. Phase power winding current 


Figurel3. Phase control winding current 
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Figure 14. DC voltage 


Figure 15. Blade pitch angle 
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Figure 16. Power coefficient Cp variation 


Figure 17. Power coefficient Cp variation 
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Figure 18. Zoom of phase PW current and 
voltage 


Figure 19. Power coefficient Cp variation 
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7. CONCLUSION 

In this study we presented in detail sensorless control strategy for (BDFG) in variable speed wind 
turbine generators used Luenberger observer, a vector control strategy using power winding flux-oriented 
scheme is proposed to assess the decoupage of active and the reactive power, the observer gains are selected 
by the pole placement method and the stability of the observer is analyzed using the Lyapunov theory. 

The simulation results show effectiveness of the optimal power sensorless operating methods in low 
and high wind speed, we can conclude the MPPT senseless operating methods proposed only by measuring 
phase voltages and currents therefore it can improve the control system dependability and energy conversion 
competence efficiency. 


Appendix 


Table 1 The Electrical Parameter of BDFG 



PW 

cw 

Rotor 

Resistance (^) 

Rp =1.732 

R ° =1.079 

Rf = 0.473 

self-inductance (mH) 

Lp =714.8 

^ =121.7 

Lr =132.6 

Mutual inductance (mH) 

Mp =242.1 

Mc =59.8 
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